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Abstract 

We evaluate the expected level of foreground contamination to the 
cosmic microwave background (CMB) polarised radiation, focusing on 
the diffuse emission from our own Galaxy. In particular, we perform a 
first attempt to simulate an all sky template of polarised emission from 
thermal dust. This study indicates that the foreground contamination to 
CMB B modes is likely to be relevant on all frequencies, and even at high 
Galactic latitudes. 

We review the recent developments in the design of data analysis tech- 
niques dedicated to the separation of CMB and foreground emissions 
in multi-frequency observations, exploiting their statistical independence. 
We argue that the high quality and detail of the present CMB observations 
represent an almost ideal statistical dataset where these algorithms can 
operate with excellent performance. We explicitly show that the recovery 
of CMB B modes is possible even if they are well below the foreground 
level, working at the arcmin. resolution at an almost null computational 
cost. This capability well represents the great potentiality of these new 
data analysis techniques, which should be seriously taken into account for 
implementation in present and future CMB observations. 



1 Introduction 

Among all cosmological observables, the Cosmic Microwave Background (CMB) 
polarisation stores cosmological perturbations in the smartest way. The Stokes 
parameters Q and U are non-locally combined to get the E and B modes, fea- 
turing opposite parity relations; as for total intensity (T), E modes are activated 
by all cosmological perturbations, and indeed the two are strongly correlated 
{TE); on the other hand, B modes are excited by vectors and tensors only; since 
no significant power is experted from vectors, a detection of B modes would be 
a strong indication of the presence of cosmological gravitational waves. For a 
comprehensive review of the CMB physics, see Hu et al. (1998) and references 
therein. 

On the observation side, we are right now in the epoch in which the CMB 
is revealing its finest structure. The Wilkinson Microwave Anisotropy Probe 
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Figure 1: Upper panels: simulated templates for the synchrotron Q Stokes 
parameter. Lower left panel: synchrotron spectral index template. Lower right 
panel: Q simulated template for the dust. 

(WMAP) satellite 1 released recently the first year CMB data, for T and TE, 
unveiling all sky anisotropics with about 20 arcmin. re solution on a frequency 
range between 22 and 90 GHz fsee lBennett et alll2003t and references therein); 
balloon-borne and ground-based T observations cover angular scales from a few 
arc minutes to tens of de grees, and also aim at measuring CMB E and B modes 
fsee iKovac et all [20021 for a first E mode detection). The Planck satellite 2 
will provide total intensity and polarisation full sky maps with 5' resolution and 
a sensitivity of a few (iK, on nine frequency channels between 30 and 857 GHz. 

This work focuses on the foreground contamination to the CMB polarised 
emission. This issue is gathering greater and greater atte ntion as the imp or- 
tance of the CMB polarisation measurements increases fsee lDe Zottil 120021 for 
reviews) . In Section[21we estimate the expected strength of foreground polarised 
emission with respect to the CMB; in Section[31we describe the recent progresses 
in the design of data analysis tools designed to reduce such a contamination; in 
Section 0] we make some concluding remarks. 

1 map.gsfc .nasa.gov 

2 astro/estec . esa.nl/SA-general/Proj ects/Planck 
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2 Foreground contamination to the CMB 



We focus on the diffuse emission from our own Galaxy. Indeed, at microwave 
frequencies, the contaminati on from Galactic, extra-Galactic radio and infrared 
sources (see iDe ZottiL 12002. and references th erein), as well as the polarised 
Sunyaev-Zel Dovich emission fsee IWhitel . I2003L proceedings of this conference), 
are likely to be relevant for the damping tail of the CMB power spectrum, i.e. at 
arcmin. angular scales, corresponding to multipolcs I of the order 1000, which 
sets the limit of validity for our analysis. 

The main diffuse foreground at low microwave frequencies, say ^ 100 GHz 
is expected from Galactic synchrotron. The available high resolution data in 
the radio band, ap proximately 10 arcmin., are at low and m edium Galactic 
latitudes, |6| < 20° founcan et all Il999t lUvaniker et al"lll999|k degree resolu- 
tion data l|Brouw fe Spoelstral . ll976(l cover about half of the sky and reach high 
latitudes. These templates have been used to evaluate the a ngular power spec- 
trum of synchrotron in the radio band, up to arcmin. scale llTucci et all l2000t 

Baccigaluni et~aTlEioillc^"rdino et allEooilTucci et alll2002HHe Oliveira-Costa et all 



proce e dings of this conference). 
Tucci et al.l l|2000l) estimated Cf ' B oc 2 , regardless of the Galactic re- 



gion considered' note t hat the same behavior w as quoted for the first time by 
iTesrmark et al.l l)200fll . iBaccigalupi et al 1 i|200ll) found the same result, regard- 
less also of the Galactic latitude, up to the highest probed \b\ ~ 20°; they also 
estimated the level of Faraday depolarisation to be important, but not such 
to mask substantially the true synchrotron emission; moreover, they estimated 
the powe r on super-degree angular scales, corresponding to multipoles I < 200, 
using thelBrouw fc Spoelstral l|l97fil) data, finding a steeper behavior, Q oc £~ 3 . 
iGiardino et alJ l|2002l) estimated the power spectrum of the cosine and sine of 
the polarisation angle up to the arcmin. scale, and out of that constructed a 
spatial template assuming Gaussian distribution; the resulting pattern was then 
used to simulate templates for Q and U out of the total intensity Haslam map 
at 408 MHz, assumed to be synchrotron dominated and theoretically polarised 
at 75%; they also estimated t he synchrotron spec tral index as inferred by multi- 
frequency radio observations. iTucci et alJ ll2002l) extended the power sp e ctrum 
analysis up to a few thousand multipoles. and Ide Oliveira-Costa et all l|2003|) 
performed a detailed study of the E and B modes in the large scale data. 

The CMB contamination from diffuse Galactic dust is much more uncer- 
tain. The main emission mechanism arises from the thermal emission of mag- 
netized dust grains, which get locally aligned by the Galactic magnetic field 
( Lazarian & Prune 1 120021) . Moreover, indications of a non-thermal dust emis- 
sion, possibly connected with rapid rotations of the grains and extending to 
much lower frequencies, has been conjectured in several w orks, although in the 
WMAP data there is no evidence of such component (see iBennett et, all I2003L 
and references therein). In total intensity, a 6' resolution all sky template of 
dust emission is available at 100 /zm, and on the basis of a correct interpreta- 
tion it can be sca l ed at microwave frequencies and compared with the CMB (see 
iFinkbeiner et allll999l and references the rein). In polarisation, a great achieve- 
ment has been obtained by Archeops (see IPonthieu et all 120031 proceedings of 
this conference), which measured a 3 — 5% polarised dust signal on the Galactic 
plane. On the basis of these two ingredients, we propose here a naively simu- 
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Figure 2: Signal distribution for the Stokes parameter Q (left) and Galactic 
latitude dependence (right) for Sq (top), Sb (middle) and dust (bottom), at 
100 GHz. 

lated template of th e polarised dust emission, which agrees remarka bly well with 
previous estimates l|Lazarian fc Prune E 120021: IPrunet et all l2000ft . We assume 
perfect correlation between total and polarised intensity and a 5% polarisation 
degree. T hen we assume the sam e polarisation angle pattern of the synchrotron 
emission l|Giardino et all 120021) . That is certainly quite drastic, but it cor- 
responds to the assumption that the Galactic magnetic field is 100% efficient 
imprinting the polarisation direction to the dust a nd synchrotron emissi on. 

The Q Stokes parameter template 3 obtained bv lGiardino et ail l|2002h . here- 
after Sg, is shown in figure H together with the spectral index template, upper 
and lower left panels, respectively; we also show the S r template obtained rescal - 
ing Sg to match the power spectrum predicted bv iBacciealupi et all l)200l|) . 

3 Unless otherwise specified, the maps and power spectra we show are in Kelvin, antenna 
temperature. 
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upper right panel; we also show the Q template for the polarised thermal dust 
emission as described above, in the right lower panel. Sg, Sb and dust are 
shown in a non-linear scale in order to highlight the behavior at high Galac- 
tic latitudes; the corresponding template for U is qualitatively equivalent; they 
have been extrapolated to 100 GHz, in antenna temper ature, using the spectra l 
index template for synchrotron and the latest recipe bv lFinkbeiner et alJ l|l999|) 
for dust. By visual inspection, the Sb model appear to be less structured on 
small angular scales with respect to the So and the dust template; the latter 
has the sharpest decrease with Galactic latitude, as we discuss below. 

Before quantifying the contamination to the CMB polarisation angular power 
spectrum, it is useful to give a look at the whole foreground statistics and sky 
distribution, plotted in figure |2 from top to bottom, the left and right panels 
represent signal distribution and root mean square as a function of b for Sg, 
Sb and dust. The signal distributions are far from Gaussianity, and exhibit 
a marked super-Gaussian behavior, as the total intensity maps from which all 
these templates were derived. Also, as we anticipated, the dust and Sb tem- 
plates present the strongest and weakest dependence on the Galactic latitude, 
respectively. For the dust, the reason is the assumed perfect correlation between 
total and polarised intensity. For Sb, as we see now more in detail, most of the 
power lies on the small angular scales, which propagate well outside the Galactic 
plane. 

In Fig. |3| we report the power spectra of the Sg, Sb and dust templates, 
compared to the CMB 4 , on all sky and with a Galactic cut excluding the region 
with latitude \b\ < 20°. The plots are at 100 GHz, in antenna temperature. The 
synchrotron frequency scaling according to the template in figure ^ is approxi- 
mately [100/f (GHz)] 5 ' 5 , making the contam ination rapidly worse at lower fre- 
quencies. For the dust, the scaling is given bv lFinkbeiner et"aTI l|l999F ). which in 
the frequencies of interest here, say the interval covered by WMAP and Planck, 
is approximately oc zA 7 /[exp [hv jkBTn) — 1], where h and ks are the Planck 
and Boltzmann constant, and the thermodynamical dust temperature Tn ~ 18 
Kelvin; the departure from a pure black bod y behavior is due to t he heating that 
the interstellar medium gets from stars fsee lPnmet et all 12000. and references 
therein) . The CMB fluctuations in antenna temperature scale as the derivative 
of the black body brightness oc v 2 exp (hu /ksTcMB) /[exp (hv/kBTcMB) — I] 2 - 
As a general feature, it can be noted how the Galactic emission has almost 
equal power on E and B modes, according to the most natural expectation for 
a non-cosmological signal (Zaldarriaga 2001). 

On a full sky analysis, the polarised dust dominates. This is likely to be the 
case if the observed 3 — 5% polarisation degree really holds for all scales and 
the polarisation angle statistics is not too far from the synchrotron one. On 
the other hand, this signal does show the most abrupt drop when the Galaxy is 
cut out, which cleans almost completely the CMB E and TE modes, remaining 
however substantial for the B ones; note th at this case agrees remarkably well 
with the level quoted bv lPrunet et"afl §000) at intermediate Galactic latitudes. 
On all sky, the synchrotron contamination is severe on large angular scales, al- 
most covering the CMB reionization bump for E and B, being somewhat less 
important for TE. On E, the Sb signal drops rapidly below the CMB at a few 

4 We realize a CMB template according to the WMAP data: h = 0.72 , Q A = 0.7 , fl b h 2 = 
0.0228 , r = 0.117 , ng = 0. 96, also assuming a te nsor to scalar ratio r = 30%, rather close to 
the current upper limit fsee lBennett et all 2003, and references therein). 



5 



tens in £, while Sq affects the first CMB E acoustic oscillation. An important 
point is that the Galactic cut is not as effective as for the dust: the reason is the 
weaker Galactic latitude dependence for synchrotron, as it is evident in figure 
El As a result, the CMB B modes are likely to be severely contaminated even 
when the Galaxy is cut out. 

We conclude that our present knowledge of the diffuse Galactic polarisation 
emission indicate that CMB TE and E modes are relatively clean, at least 
at medium and high Galactic latitudes and at frequencies around 100 GHz, 
while the B modes are likely to be severely contaminated, at all frequencies and 
Galactic latitudes. 

3 Component Separation 

Despite of the severe foreground contamination, the CMB B modes are expected 
to carry most important details on the nature of cosmological perturbations, as 
we stressed in the Introduction. For this reason it is worth to study data analysis 
techniques capable to clean the CMB from foreground contamination. This is 
feasible by exploiting the statistical differences of CMB and diffuse foregrounds, 
starting from a wide enough multi-frequency coverage of a given observation. 
The data can be sketched as 

x(rif,n p ) = A(n/, n c ) ■ s(n c , n p ) + n(ra/, n p ) , (1) 

where x represents the data at the nj frequencies, s the n c components cor- 
responding to each emission, and n the instrumental noise; A is the mixing 
matrix, scaling and mixing the emissions s at any observed frequency, while 
n p is the number of pixels of the observation, or equivalently the number of 
harmonic coefficients if the analysis is performed in the spectral domain. Note 
that in general the angular resolution varies among the frequency channels. The 
target is to recover s from x; as we saw in the previous Section, the problem 
is undetermined since we know A only partially, and the guesses on s are not 
reliable for stabilizing the inversion of the relation , either if the noise statis- 
tical properties are very well known, or even in a noiseless case. On the other 
hand, a well established expectation is that the CMB and foreground emission 
are statistically independent. It has been recently shown that this feature can 
be exploited to design "blind" component separation tools, i.e. able to recover 
s even without any prior on s or A; a first algor ithm of this class was pr oposed 
for use in astrophysical component separation bv lBaccigalupi e~ ]d](|2000) based 
on the novel concept in signal processing science, the Independent Component 
Analysis: that works by finding a W(n c , rif) matrix which applied to x(n/, n c ) 
in (JJJ brings to s + Wn. The criterion to find W is to look for the statistically 
independent patterns in x: that is clearly not unique and can be specialized for 
the problem at hand. A fa st version (FastIC A) adapted for satellite observa- 
tions has been proposed bv lMaino et alJ ((2002) ; this technique has been further 
developed to account for any av ailable and reliable prior, and successfully ap- 
plied to the real COBE data bv lMaino etafl <^2003^ . Recently, a different and 
promising blin d technique ba sed on a spectral matching of real and modeled 
data lDelabrouille et al.l l|200. c J l. 

The first applic ation of componen t separ ation to polarisation radiation has 
been performed bv iBaccigaluoi et alJ l|2003j) . with a FastICA approach: sep- 
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Figure 3: Expected angular power spectrum from Sq (left panels), Sb (middle 
panels) and dust (right panels) with respect to the CMB (dotted); in all panels, 
foreground contamination is calculated on all sky and with a Galactic cut for 
|6| < 20°. 
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aration is conducted on Q and U Stokes parameters separately or considered 
as a single signal provided that they obey the same statistics; we refer to that 
work for any detail on the algorithm procedure. Here we show a very simple 
application, with rather impressive results, to highlight the potentialities of this 
technique. We consider a noiseless mixture of polarised CMB and synchrotron 
Sq, scaling uniformly in frequency, at 50 and 80 GHz, on all sky and pixel size of 
about 3.5 arcmin.; note that we do not refer to any operating or planned exper- 
iment, and choose the frequencies among those in which CMB and synchrotron 
are roughly comparable. The two channels are input to the FastICA to re- 
cover synchrotron and CMB separately. Figure 0] show the input sky angular 
power spectra for E (left) and B (right) together with the spectra of the recov- 
ered CMB as output by FastICA. As it is evident, a CMB major cleaning was 
achieved: the dotted line is the input sky B power, while the recovered CMB 
(solid) line is actually almost indistinguishable from the input one (dashed). We 
stress that this is without the use of any prior, i.e. not knowing any piece of A 
or s. The algorithm knows only the data at the two frequencies, and looks for 
the most statistically independent components in those data. In addition, ap- 
plying this technique in particular is computationally free: the results in figure 
were obtained in about five minutes on an ordinary single processor personal 
computer. 

In particular, the CMB B modes have been impressively cleaned: interpreting 
the foreground emission as noise, the signal to noise ratio is 0.1 or less, at all 
multipoles, as it can bee seen in the right panel of the figure. CMB E modes 
have also been cleaned wherever the foreground contamination was relevant. 
The algorithm outputs also A, which differs from the real one by percent or 
less. 

The reason of such an excellent performance is just the high level of detail of 
the underlying data. At arcmin. resolution, the maps are almost ideal datasets 
where the algorithm can look for the least statistically linked templates present 
into the maps themselves. Even if all the instrumental details were neglected, 
what we show here should be enough to gather attention on these new data 
analysis techniques, as we further stress in the next Section. 

4 Concluding remarks 

In this work we show that at frequencies around 100 GHz the diffuse Galactic 
polarisation emission is likely to be rather weak in comparison to the CMB 
TE and E modes, at least at medium and high Galactic latitudes, while the 
B ones are seriously contaminated on all sky. At lower and higher frequencies, 
the contamination gets rapidly worse because of the steep frequency behavior 
of synchrotron and dust, respectively. 

We stress that the new concepts in signal processing science can be exploited 
in this context, to cancel or greatly reduce such a contamination. In particu- 
lar, the statistical independence of CMB and foreground emission is suitable 
to achieve this goal: we give an impressive exam ple based on the FastICA 
technique in polarisation l|BaccigaluDi et all [2003). where CMB B modes are 
recovered out of a foreground emission about ten times stronger, without any 
prior information but their statistical independence, and at an almost null com- 
putational cost. 
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Figure 4: E and B modes of the original sky at 50 and 80 GHz (respectively high 
and low dotted curves) and recovered CMB E and B modes after component 
separation (solid), at 80 GHz; the true CMB (dashed), almost indistinguishable 
from the recovered one, is also plotted. 



These features represent the "nominal" capabilities of these algorithms, i.e. 
their performance on ideal datasets: cleaning the cosmological signal from sub- 
stantial foreground contamination, exploiting only their natural diversity. This 
is in our opinion a great achievement, providing the basis on which the im- 
plementation to any particular observation can be built. The latter step, of 
course, presents challenging issues: the noise level, sky distribution and non- 
stationarity, as well as the non-uniform spectral index of the components to 
recover, are good examples. However, we have already the indications of a sub- 
stantial stabil ity of the nominal perf ormance against instrume ntal features and 
systematics: iMaino et all l)2002j) and iBaccieraluni et al.l I^OO.^ . respectively for 
total intensity and polarisation, perform a whitening of simulated Planck data 
in order to get rid of the n oise bulk in the Fa stICA separation process. Most 
importantly, the results of IMaino et alJ l|200. c {f > show that these techniques are 
flexible enough to be applied to real data, the COBE ones, getting outstand- 
ing results, namely the correct large scale CMB power spectrum amplitude and 
shape. 

While waiting for the era of large scale CMB polarisation measurements, the 
next obvious target is a successful application to the WMAP data, meaning the 
recovery of CMB data on sky regions affected by Galactic contamination which 
are normally cut out. 

The author warmly thanks Gianfranco De Zotti, Giovanna Giardino, Eric 
Hivon, Francesca Perrotta, George Smoot and Radek Stompor for useful discus- 
sions. The HEALPix pixelisation scheme, available at www.eso.org/healpix, 
by A.J. Banday, M. Bartelmann, K.M. Gorski, F.K. Hansen, E.F. Hivon, and 
B.D. Wandelt, has been extensively used. 
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